Discovery of a large set of SNP and SSR genetic markers by high-throughput sequencing of pepper (Capsicum annuum)

Genetic markers based on single nucleotide polymorphisms (SNPs) are in increasing demand for genome mapping and fingerprinting of breeding populations in crop plants. Recent advances in high-throughput sequencing provide the opportunity for whole-genome resequencing and identification of allelic variants by mapping the reads to a reference genome. However, for many species, such as pepper (Capsicum annuum), a reference genome sequence is not yet available. To this end, we sequenced the C. annuum cv. "Yolo Wonder" transcriptome using Roche 454 pyrosequencing and assembled de novo 23,748 isotigs and 60,370 singletons. Mapping of 10,886,425 reads obtained by the Illumina GA II sequencing of C. annuum cv. "Criollo de Morclos 334" to the "Yolo Wonder" transcriptome allowed for SNP identification. By setting a threshold value that allows selecting reliable SNPs with minimal loss of information, 11,849 reliable SNPs spread across 5919 isotigs were identified. In addition, 853 single sequence repeats were obtained. This information has been made available online

Analysis of the genetic diversity of influenza A viruses using next-generation DNA sequencing

Background: Influenza viruses exist as a large group of closely related viral genomes, also called quasispecies. The composition of this influenza viral quasispecies can be determined by an accurate and sensitive sequencing technique and data analysis pipeline. We compared the suitability of two benchtop next-generation sequencers for whole genome influenza A quasispecies analysis: the Illumina MiSeq sequencing-by-synthesis and the Ion Torrent PGM semiconductor sequencing technique. Results: We first compared the accuracy and sensitivity of both sequencers using plasmid DNA and different ratios of wild type and mutant plasmid. Illumina MiSeq sequencing reads were one and a half times more accurate than those of the Ion Torrent PGM. The majority of sequencing errors were substitutions on the Illumina MiSeq and insertions and deletions, mostly in homopolymer regions, on the Ion Torrent PGM. To evaluate the suitability of the two techniques for determining the genome diversity of influenza A virus, we generated plasmid-derived PR8 virus and grew this virus in vitro. We also optimized an RT-PCR protocol to obtain uniform coverage of all eight genomic RNA segments. The sequencing reads obtained with both sequencers could successfully be assembled de novo into the segmented influenza virus genome. After mapping of the reads to the reference genome, we found that the detection limit for reliable recognition of variants in the viral genome required a frequency of 0.5% or higher. This threshold exceeds the background error rate resulting from the RT-PCR reaction and the sequencing method. Most of the variants in the PR8 virus genome were present in hemagglutinin, and these mutations were detected by both sequencers. Conclusions: Our approach underlines the power and limitations of two commonly used next-generation sequencers for the analysis of influenza virus gene diversity. We conclude that the Illumina MiSeq platform is better suited for detecting variant sequences whereas the Ion Torrent PGM platform has a shorter turnaround time. The data analysis pipeline that we propose here will also help to standardize variant calling in small RNA genomes based on next-generation sequencing data

Characterization of genome-wide ordered sequence-tagged Mycobacterium mutant libraries by Cartesian Pooling-Coordinate Sequencing

Reverse genetics research approaches require the availability of methods to rapidly generate specific mutants. Alternatively, where these methods are lacking, the construction of pre-characterized libraries of mutants can be extremely valuable. However, this can be complex, expensive and time consuming. Here, we describe a robust, easy to implement parallel sequencing-based method (Cartesian Pooling-Coordinate Sequencing or CP-CSeq) that reports both on the identity as well as on the location of sequence-tagged biological entities in well-plate archived clone collections. We demonstrate this approach using a transposon insertion mutant library of the Mycobacterium bovis BCG vaccine strain, providing the largest resource of mutants in any strain of the M. tuberculosis complex. The method is applicable to any entity for which sequence-tagged identification is possible

WAVELET BASED FUNCTIONAL MODELS FOR TRANSCRIPTOME ANALYSIS WITH TILING ARRAYS

For a better understanding of the biology of an organism a complete description is needed of all regions of the genome that are actively transcribed. Tiling arrays can be used for this purpose. Such arrays allow the discovery of novel transcripts and the assessment of differential expression between two or more experimental conditions such as genotype, treatment, tissue, etc. Much of the initial methodological efforts were designed for transcript discovery, while more recent developments also focus on differential expression. To our knowledge no methods for tiling arrays are described in the literature that can both assess transcript discovery and identify differentially expressed transcripts, simultaneously. The wavelet based functional model developed in this paper is designed to fill this methodological void. As opposed to existing methods, our statistical framework also permits a natural integration of preprocessing into the standard statistical analysis flow of tiling array data. We use Johnson transformations, which are based on cumulants, for computing false discovery rates (FDRs) and Bayesian credibility intervals for the estimates of the effect functions within the data space. A case study illustrates that our model is well suited for a simultaneous assessment of transcript discovery and differential expression, while remaining competitive with methods that perform only one of these tasks

LPS resistance of SPRET/Ei mice is mediated by Gilz, encoded by the Tsc22d3 gene on the X chromosome

Natural variation for LPS-induced lethal inflammation in mice is useful for identifying new genes that regulate sepsis, which could form the basis for novel therapies for systemic inflammation in humans. Here we report that LPS resistance of the inbred mouse strain SPRET/Ei, previously reported to depend on the glucocorticoid receptor (GR), maps to the distal region of the X-chromosome. The GR-inducible gene Tsc22d3, encoding the protein Gilz and located in the critical region on the X-chromosome, showed a higher expressed SPRET/Ei allele, regulated in cis. Higher Gilz levels were causally related to reduced inflammation, as shown with knockdown and overexpression studies in macrophages. Transient overexpression of Gilz by hydrodynamic plasmid injection confirmed that Gilz protects mice against endotoxemia Our data strongly suggest that Gilz is responsible for the LPS resistance of SPRET/Ei mice and that it could become a treatment option for sepsis

Physiological and transcriptomic evidence for a close coupling between chloroplast ontogeny and cell cycle progression in the pennate diatom Seminavis robusta

Despite the growing interest in diatom genomics, detailed time series of gene expression in relation to key cellular processes are still lacking. Here, we investigated the relationships between the cell cycle and chloroplast development in the pennate diatom Seminavis robusta. This diatom possesses two chloroplasts with a well-orchestrated developmental cycle, common to many pennate diatoms. By assessing the effects of induced cell cycle arrest with microscopy and flow cytometry, we found that division and reorganization of the chloroplasts are initiated only after S-phase progression. Next, we quantified the expression of the S. robusta FtsZ homolog to address the division status of chloroplasts during synchronized growth and monitored microscopically their dynamics in relation to nuclear division and silicon deposition. We show that chloroplasts divide and relocate during the S/G2 phase, after which a girdle band is deposited to accommodate cell growth. Synchronized cultures of two genotypes were subsequently used for a cDNA-amplified fragment length polymorphism-based genome-wide transcript profiling, in which 917 reproducibly modulated transcripts were identified. We observed that genes involved in pigment biosynthesis and coding for light-harvesting proteins were up-regulated during G2/M phase and cell separation. Light and cell cycle progression were both found to affect fucoxanthin-chlorophyll a/c-binding protein expression and accumulation of fucoxanthin cell content. Because chloroplasts elongate at the stage of cytokinesis, cell cycle-modulated photosynthetic gene expression and synthesis of pigments in concert with cell division might balance chloroplast growth, which confirms that chloroplast biogenesis in S. robusta is tightly regulated

Molecular phenotyping of the pal1 and pal2 mutants of Arabidopsis thaliana reveals far-reaching consequences on phenylpropanoid, amino acid, and carbohydrate metabolism

The first enzyme of the phenylpropanoid pathway, Phe ammonia-lyase (PAL), is encoded by four genes in Arabidopsis thaliana. Whereas PAL function is well established in various plants, an insight into the functional significance of individual gene family members is lacking. We show that in the absence of clear phenotypic alterations in the Arabidopsis pall and pal2 single mutants and with limited phenotypic alterations in the pall pal2 double mutant, significant modifications occur in the transcriptome and metabolome of the pal mutants. The disruption of PAL led to transcriptomic adaptation of components of the phenylpropanoid biosynthesis, carbohydrate metabolism, and amino acid metabolism, revealing complex interactions at the level of gene expression between these pathways. Corresponding biochemical changes included a decrease in the three major flavonol glycosides, glycosylated vanillic acid, scopolin, and two novel feruloyl malates coupled to coniferyl alcohol. Moreover, Phe overaccumulated in the double mutant, and the levels of many other amino acids were significantly imbalanced. The lignin content was significantly reduced, and the syringyl/guaiacyl ratio of lignin monomers had increased. Together, from the molecular phenotype, common and specific functions of PAL1 and PAL2 are delineated, and PAL1 is qualified as being more important for the generation of phenylpropanoids

Genome-wide screening for cis-regulatory variation using a classical diallel crossing scheme

Large-scale screening studies carried out to date for genetic variants that affect gene regulation are generally limited to descriptions of differences in allele-specific expression (ASE) detected in vivo. Allele-specific differences in gene expression provide evidence for a model whereby cis-acting genetic variation results in differential expression between alleles. Such gene surveys for regulatory variation are a first step in identifying the specific nucleotide changes that govern gene expression differences, but they leave the underlying mechanisms unexplored. Here, we propose a quantitative genetics approach to perform a genome-wide analysis of ASE differences (GASED). The GASED approach is based on a diallel design that is often used in plant breeding programs to estimate general combining abilities (GCA) of specific inbred lines and to identify high-yielding hybrid combinations of parents based on their specific combining abilities (SCAs). In a context of gene expression, the values of GCA and SCA parameters allow cis- and trans-regulatory changes to be distinguished and imbalances in gene expression to be ascribed to cis-regulatory variation. With this approach, a total of 715 genes could be identified that are likely to carry allelic polymorphisms responsible for at least a 1.5-fold allelic expression difference in a total of 10 diploid Arabidopsis thaliana hybrids. The major strength of the GASED approach, compared to other ASE detection methods, is that it is not restricted to genes with allelic transcript variants. Although a false-positive rate of 9/41 was observed, the GASED approach is a valuable pre-screening method that can accelerate systematic surveys of naturally occurring cis-regulatory variation among inbred lines for laboratory species, such as Arabidopsis, mouse, rat and fruitfly, and economically important crop species, such as corn

Sequence-specific protein aggregation generates defined protein knockdowns in plants

Protein aggregation is determined by short (5-15 amino acids) aggregation-prone regions (APRs) of the polypeptide sequence that self-associate in a specific manner to form beta-structured inclusions. Here, we demonstrate that the sequence specificity of APRs can be exploited to selectively knock down proteins with different localization and function in plants. Synthetic aggregation-prone peptides derived from the APRs of either the negative regulators of the brassinosteroid (BR) signaling, the glycogen synthase kinase 3/Arabidopsis SHAGGY-like kinases (GSK3/ASKs), or the starch-degrading enzyme alpha-glucan water dikinase were designed. Stable expression of the APRs in Arabidopsis (Arabidopsis thaliana) and maize (Zea mays) induced aggregation of the target proteins, giving rise to plants displaying constitutive BR responses and increased starch content, respectively. Overall, we show that the sequence specificity of APRs can be harnessed to generate aggregation-associated phenotypes in a targeted manner in different subcellular compartments. This study points toward the potential application of induced targeted aggregation as a useful tool to knock down protein functions in plants and, especially, to generate beneficial traits in crops

A screening assay for Selective Dimerizing Glucocorticoid Receptor Agonists and Modulators (SEDIGRAM) that are effective against acute inflammation

It has been suggested that glucocorticoid receptor (GR) agonists that promote GR homodimerization more than standard glucocorticoids such as Dexamethasone could be more effective anti-inflammatory molecules against acute and life-threatening inflammatory conditions. To test this hypothesis, we set up a screening pipeline aimed at discovering such Selective Dimerizing GR Agonists and Modulators (SEDIGRAM). The pipeline consists of a reporter gene assay based on a palindromic glucocorticoid responsive element (GRE). This assay represents GR dimerization in human A549 lung epithelial cells. In the pipeline, this is followed by analysis of endogenous GRE-driven gene expression, a FRET assay confirming dimerization, and monitoring of in vitro and in vivo anti-inflammatory activity. In a proof of principle experiment, starting from seven candidate compounds, we identified two potentially interesting compounds (Cortivazol and AZD2906) that confer strong protection in a mouse model of aggressive TNF-induced lethal inflammation. A screening pipeline for SEDIGRAM may assist the search for compounds that promote GR dimerization and limit overwhelming acute inflammatory responses